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EFFECT OF GAGE SIZE ON THE MEASUREMENT OF LOCAL HEAT FLUX 
by Kenneth J. Baumeister and S. Stephen Papell 
Lewis Research Center 


SUMMARY 

General formulas are derived for determining gage averaging errors of strip-type 
heat flux meters used in the measurement of one- dimensional heat flux distributions. 

The local averaging error e(x) is defined as the difference between the measured value 
of the heat flux and the local value which occurs at the center of the gage. The present 
report will develop some simple expressions for estimating e(x) by the use of a trun- 
cated Taylor series approximation. In the analysis, the gage is assumed to operate 
under steady-state conditions and to be "perfect"; that is, it is assumed to measure 
exactly the heat flux that falls on it. In actual practice, however, conduction and radia- 
tion losses must be accounted for by calculation or calibration. In addition, a correction 
procedure is presented which allows a better estimate for the true value of the local heat 
flux. As an example of the technique, the formulas are applied to the cases of heat 
transfer to air slot jets impinging on flat and concave surfaces, as might occur in the 
internal cooling channels of turbine blades. One -dimensional heat flux distributions oc- 
cur when cooling is accomplished with slot jets or with round hole jets aligned in tightly 
packed rows. It has been found that for many practical problems, the use of very small 
heat flux gages is often unnecessary. 


INTRODUCTION 

The aviation and power-generation industries have been concerned with heat flux 
measurements (refs. 1 and 2). Prior to 1960, measurements of heat flux (ref. 3) were 
generally concerned with average values taken over large areas. However, because of 
interest in prediction of local temperature gradients and thermal stresses, such as in 
the thermal design of turbine blades, recent work has been concerned with local heat 
flux measurements. 

In choosing a particular heat flux meter, such as a slug-type, Gardon, or ablation 
calorimeter, the experimenter is faced with the choice of gage size. The gage can be 



assumed to measure the average heat flux across its face, rather than a local value at 
its center (see fig. 1). In principle, the smallest possible heat flux gage should be used 
for accurate local heat flux measurements. However, the smaller the size of the gage, 
the greater will be the fabrication, handling, and calibration problems, as well as cost. 
This is illustrated in figure 2. The averaging error e(x) (difference between center 
value and average measured value of the heat flux) will decrease for smaller gage sizes, 
as illustrated in figure 2. On the other hand, for very small gages, calibration errors 
may increase and costs will rise because of fabrication difficulties. Clearly, some opti- 
mum choice of gage will generally be required, depending on the heat flux distribution 
and fabrication and calibration difficulty. The role of the heat flux distribution will be 
discussed in detail in the main body of this report. 

To estimate the optimum size heat flux gage, the experimenter would like to know 
the "averaging" error e(x). At the present time, however, the literature is void of any 
means of predicting e(x) for these gages. In the present report, some relatively simple 
expressions will be developed for estimating e(x) by the use of a truncated Taylor series 
approximation. With the resulting analytical expressions, the experimenter can conve- 
niently estimate the maximum gage size permissible for a given acceptable experimental 
error, provided he knows the general shape of the heat flux distribution curve. In detail, 
the purpose of this report is as follows: 

(1) To present an expression for the error e(x) due to gage averaging for one- 

dimensional heat flux distributions as measured with a strip-type gage (The 
gage is assumed to operate under steady- state conditions and to be "perfect"; 
i. e. , it is assumed to measure exactly the heat flux that falls on it. In actual 
practice, however, there are losses, such as thermal conduction leakage or 
thermal radiation losses. These must be accounted for by calculation or cali- 
bration. ) 

(2) To present parametric plots of the error distribution e(x) for some common 

types of heat flux distributions that might occur in practice 

(3) To apply the analysis to some problems in jet impingement heat transfer 

(4) To present a means of correcting measured heat flux data to account for gage 

size 

SYMBOLS 


a 

constant 

B 

width of slot jet (ref. 4) 

C , Cq, C p C2 

constants 

D 

diameter of hole 
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e probe averaging error 

H distance between holes 

h heat transfer coefficient 

L characteristic length (see eq. (2)) 

Nu Nusselt number 

Pr Prandtl number 

q heat flux 

<5 measured "averaged" heat flux 

in 6db 

Re Reynolds number 

AT difference between wall temperature and bulk temperature 

x distance from stagnation point 

Ax probe size (see fig. 5) 

Z n distance from slot jet to surface (ref. 4) 

j3 constant (see eq. (20)) 

rj dimensionless distance, x/L 

A 77 dimensionless width of probe, Ax/L 


ERROR ANALYSIS 
Background 

Some impinging jet experiments for rows of holes and slot jets will now be discussed 
to illustrate some of the problems involved with the measurement of the local heat flux. 
Also, as an illustrative example, the analytical techniques will be applied to some jet 
impingement data. 

Gardon and Akfirat (ref. 4) measured the local heat transfer coefficients under im- 
pinging axisymmetrical and slotted jets. Figure 3 shows some of their typical experi- 
mental results. They used a small 0.9-millimeter-diameter Gardon-type heat flux 
meter. Because of the small size of their meter, we would naturally expect very good 
local measurements. On the other hand, Tabakoff and Clevenger (ref. 5) used relatively 
large (3.8-cm wide) heat flux measurement devices (electrical heaters) to measure their 
heat transfer coefficient distribution. The validity of their local measurements is cer- 
tainly subject to question. Figure 4 shows atypical experimental result from reference 5 
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for rows of holes impinging on a concave surface. The question arises, how much de- 
tailed local information can we deduce from the measured heat transfer coefficient dis- 
tribution? Of course, one obvious problem with such gages, as used in reference 5 , is 
the possibility that finely structured peaks might occur and be undetected. Obviously, 
such measurements would lead to greatly underestimated temperature and thermal stress 
distributions if peaks occurred. 

The manner in which the gage is used is also extremely important in sizing the gage. 
If the gage can be positioned at all values of x, as was done by Gardon and Akfirat in 
reference 4 by moving their plate and heat flux gage relative to a fixed jet, the local 
peaks have a much better chance of being observed. On the other hand, if for a particu- 
lar experiment the gages must be placed in a fixed position, as was done in reference 5 
for the data shown in figure 4 herein, then it will not be possible to determine if any 
local peaks occur within one gage size. Also, for very small fixed gages, local peaks 
between the gages will also go undetected. A local peak, for example, might occur if the 
boundary layer changes from laminar to turbulent flow, that is, if transition occurs. 

In addition to sizing gages for new experiments, there is also a need to judge past 
work. For example, how accurate was Gardon and Akfirat's measurement of the peak 
heat transfer coefficient (fig. 3 (a)), or how much higher is the local heat transfer coeffi- 
cient at the stagnation point than the average value reported by Tabakoff and Clevenger 
(fig. 4 )? 


Averaging Technique 

Consider the one- dimensional heat flux distribution q(x), as shown in figure 5 , 
which is being measured by a heat flux meter of size Ax. For the purpose of this re- 
port, we assume that the gage measures exactly the total heat flux that falls on it; that 
is, we assume that we can account for all losses from the heat flux meter, such as tem- 
perature gradients along the wall, and correct for them. In this section, we are con- 
cerned only with the problem of a large probe diameter washing out local effects. 

The average heat flux measured by the gage can be expressed as 


hneas 


( V)= — f 

An •'n- 


77+ A 77/2 


A 77 •'77- A 77/2 

where 77 is a dimensionless distance chosen as 


q(77)d77 


V = —i 
L 


A77 = 


Ax 


( 1 ) 


( 2 ) 
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Thus, the error in the measurement of the local heat flux can be expressed as 


e(v) = 


^meas^ 


/ ^ At ?/2 
A tj /2 


q(r?)d?7 


<1(7?) 


Ar?q(r?) 


( 3 ) 


By using equation ( 3 ), we can evaluate the errors associated with some common q (77) 
distributions for various dimensionless probe sizes A 77. 


General Case 


Before we begin investigating specific heat flux distributions, let us consider proper- 
ties of q( 77) which could lead to large measurement errors. The function q(7?) can be 
expanded in a Taylor series about any point "a" as follows: 

q(?7) = q(a) + (77 - a)q I (a) + ^ ~ a ^ - q n (a) + q ni (a) + q IV (a) + . . . ( 4 ) 

2 ! 3 ! 4 ! 


where (I, II, . . . ) represents derivatives with respect to 77. 

Substituting equation ( 4 ) into equation ( 3 ) and performing the integration gives 


e(,) = aW-i + s I w^-a) + a^(a) 

q(77) q(rf) 2q(??) 


(77- a) 2 



+ ^ 4(77 - a)** + At7 2 (t 7 - a) 
24 q( 77 ) L 


q IV (a) 

24 


^ .At? 0 ^2 , A77 


(t? - a) + — z_ (77 - a) + 


80 


( 5 ) 


The location of the parameter "a" could be at some fixed point in the coordinate 
system, such as t? = 0 , or at the center of the heat flux gage. The simplest expression 
results when we expand about the center of the gage (a = 77). In this case, equation ( 5 ) 
becomes 


e(T7 ) _ q n (? 7 ) At ? 2 [ q IV ( 77) At ? 4 f 
q( 77) 24 q( 77) 1920 


( 6 ) 


or 


e( 7 j) = <&> ±1 

q( 77) 24 


i + ^m ±I + . . . 

q%) 80 


( 7 ) 
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For many practical problems, the contribution of the fourth and higher order terms 
in equation (5) may be neglected. That is, 


A(v) 

q H ( v) 80 


« 1 


( 8 ) 


Thus, the error in the heat flux meter can be expressed as 

e(,)»S&!AiL (9) 

q(? 7 ) 24 

The error, therefore, is proportional to the square of the size of the heat flux meter and 
to the second spatial derivative of heat flux. 

Equation (9) can now be solved to determine the size of the heat flux gage for a given 
error: 


A 1) = 




2 4 e(rj)q(ri) 

Av) 


( 10 ) 


For calculational purposes, the heat flux q(r]) is often expressed in terms of a heat 
transfer coefficient h( 77 ) as follows: 


q(??) = h( 77 )AT( 7 ?) (11) 

Substituting equation (11) into equation (9) gives the following for the error: 

+ ^iW (12) 

\ h di? a " dr < iT hi? 1 24 

The spatial temperature derivatives will, of course, depend on the particular flow situa- 
tion being measured. 


Peaks 

The expression for the error just derived is inconvenient to use when any part of the 
heat flux gage crosses a peak, such as might occur at a stagnation point, because of the 
large number of terms that would be needed to express a continuous function for the peak. 
In the previous analysis, it was assumed that the odd derivatives were continuous across 
the gage center, at 77 = 0, and thus did not contribute in the integration of equation (3). 
However, at a flux peak, symmetry is now assumed to exist for all terms in the Taylor 
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expansion. Hence, equation (3) becomes 


e(0) = 


/*A 77/2 

2 / q(?7)d77 

A 1 

At? q(0) 


(13) 


For this special case, the parameter ”a" in equation (4) takes on the value of 
zero, and the Taylor series reduces to a Maclaurin series. Substituting the Maclaurin 
series into equation (13) and performing the specified integration gives 


e (o) = q T (0) A_5 + q n (0) A?? 2 + q rn (0 ) At? 3 + g IV (0) A?? 4 + 
q(0) 4 q(0) 24 q(0) 192 q(0) 1920 


or 


e(0) = 5J0) 

q(o) 4 


q H (0) At? . q m (0) At ? 2 q IV At ? 3 


q T (0) 


q T (0) 


48 


q T (0) 480 


(14) 


(15) 


As can be seen from equation (15), both the odd and even derivatives contributed to 

the error, in contrast to equation (7), in which only the even order terms contributed to 

2 

the error. Consequently, the estimate to order At? is 

e(0) » Atj + (16 

q(0) 4 q(0) 24 


Note, if the heat flux is expanded in terms of even functions such that the first de- 
rivatives are zero at the origin, then equations (9) and (16) are identical. 

Now, we will apply the general theory to two simple example profiles. Let us as- 
sume that the heat flux distribution can be approximated by either a second-order poly- 
nomial or an exponential distribution. These two distributions were chosen because they 
can be used to approximate many practical heat flux distributions. Also, we can illus- 
trate a case where the truncated equations cannot be used, and we must rely on an exact 
integration of equation (3) or (13) to determine e(x). 


Polynomial 

An expression for the heat flux can be written in terms of a polynomial of the form 

q(7?) = C Q + C lV + C 2 t? 2 (17) 
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The preceding second-order polynomial should be sufficient for describing many mea- 
sured heat flux distributions. On the other hand, if more terms are needed in the poly- 
nomial, the extension is straightforward. 

Since the third and higher order derivatives of equation ( 17 ) are zero, the truncated 
series expressions for the error are exact. Thus, the errors are 


and 


C 9 At ? 2 /12 

e(r?) = 

C Q + C-^ 7 ] + C2V 


for 77 > -^-5 
2 


( 18 ) 


e(0) 


Ar7 + _^2_ Ar 2 

C 0 4 12 C Q 


for 77 = 0 


( 19 ) 


Exponential 

The exponential distribution for the heat flux can be expressed as 

q(7?) = Ce'^ (20) 

Substituting equation ( 20 ) into equations ( 3 ) and ( 13 ) and solving for the error distribu- 
tions gives 

e(rj) = — — sinh ^ A 77 - 1 for 77 > -^ (21) 

/3 A77 2 2 

e(0) = — — (1 - e "^ AT7/2 ) - 1 for 77=0 (22) 

j 3 A77 

If we use the truncated series expressions, equations ( 9 ) and ( 16 ), the error distribution 
becomes 


e(v) = — for 77 > ^ 
24 2 


_ -/I A77 + ff 2 At? 2 
4 24 


e(0) = M for 77 = 0 


( 23 ) 

( 24 ) 
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In equation (24), we have kept the first two terms in the series, since this will make both 
equations (23) and (24) accurate to at least order A 77 squared. 

The choice of the exponential representation is convenient, since the error is inde- 
pendent of position 77 , which was not the case for the second-order polynomial. Thus, a 
single plot will suffice for a complete mapping of the errors associated with a heat flux 
distribution, which can be approximated by an exponential. 

Parametric plots for the errors involved with the exponential approximation are dis- 
played in figures 6 and 7. Figure 6 considers measurements for 77 > A 77 / 2 , while fig- 
ure 7 considers measurements at a point of symmetry. As seen in these figures, for 
/3 values less than 2, the truncated series expression for the error can be used without 
significant error. For /3 values greater than 2, an exact integration is required. How- 
ever, for most practical problems the value of |3 is much less than 2. 

Next, we shall use the results of this section to estimate the error involved in the 
measurements of Gardon and Akfirat (ref. 4) and of Tabakoff and Clevenger (ref. 5). 


IMPINGING JET EXAMPLES 

The results of the last section will be applied to the problem of heat transfer from 
impinging jets. In particular, we will examine the data of Gardon and Akfirat (ref. 4) 
for a two-dimensional jet, and the data of Tabakoff and Clevenger (ref. 5) for a row of 
hole jets impinging on a curved surface. In both experiments, the heat flux is essenti- 
ally one- dimensional. The data of Tabakoff and Clevenger, shown in figure 3, are of 
particular interest, since we wish to examine the effect of very large heat flux meters 
on local measurements. 

Tabakoff and Clevenger used strip-type heat flux meters; consequently, the analysis 
just presented is directly applicable to their data. On the other hand, Gardon and 
Akfirat used small cylindrical meters. As a first approximation, we assume the present 
analysis can be applied to such data. 

Gardon and Akfirat's data for the slot jet (Z n /B = 8 , see fig. 3(a)) near the stagna- 
tion point will be examined first. In the vicinity of the stagnation point, the heat transfer 
coefficient can be fitted by the equation 

h( 77 ) = 85 - 13 ?? 2 for 77 ^ 1 (25) 

where 

V = ~ (26) 

B 

We will assume that the wall-to-bulk temperature gradients are negligible; thus, the 
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heat transfer coefficient and the heat flux will have a one-to-one correspondence in the 
expressions for the error, and they may be used interchangeably. Thus, if equation (25) 
is combined with equations (17) and (19), the expression for the error becomes 

e( r?) = -0.0127 A?? 2 (27) 

The 0.9-millimeter-diameter gage used by Gardon and Akfirat had a A 77 value of 
0.282. If the value of 0.282 is substituted into equation (27), the error in the absolute 
value of the heat transfer coefficient is calculated to be -0.1 percent, which is quite ac- 
ceptable. At all other points in the flow field, the error will be negligible. Had the gage 
diameter been increased by one order of magnitude, from 0.9 to 9.0 millimeter, the 
error would be about 10 percent. 

As shown in figure 8, the data (fig. 4) of Tabakoff and Clevenger (ref. 5) can be 
fitted by a polynomial of the form 

h«: 2.65 - 3.7 77 + 1.8 r? (28) 


where 


(2 x heater width) 


(29) 


Since only three measured values are presented to the left of the stagnation point, these 
values are fitted exactly by a second-order polynomial. The error at the stagnation 
point in this case becomes 

e(rj) = -0. 35 A77 + 0. 057 A?? 2 (30) 

In this particular case, A 77 equals 0.5, and the error in the absolute value at the stagna- 
tion point is -16 percent. 

The assumption was made here that the heat transfer coefficient is a monotonically 
decreasing function for 77 < 1/2. Had any singularities or peaks existed in this range, 
they would have gone undetected. 

The error at the second station in figure 4 can be evaluated from equation (18) to be 
3 percent. If the heat flux is monotonically decreasing, we see that the error in the 
measurement for this relatively large gage is quite acceptable. Also, if the gage were 
movable instead of being fixed, some peaks could even be detected by this gage. 
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CORRECTION PROCEDURE 


The application of corrections to raw data is commonplace in engineering. The tem- 
perature distributions on the insides of tubes are often determined from measurements 
taken on the outsides of tubes by appropriate consideration of the thermal conductivities 
of the tubes and other operating parameters. In a similar manner, the heat flux mea- 
surements can also be corrected to account for the error due to probe averaging. 

A better estimate for the "true" heat flux can be determined from the measured heat 
flux by solving the simple integral equation 

i /•7P-A7J/2 

— / q(? 7 )dr? = q ( 77 ) (31) 

Atj A77/2 meaS 

For the exponential and polynomial cases, this integral reduces to 

q mpas (t)) 

q(r?) = m ? aS -- (32) 

e(r?) + 1 

where the various expressions for e(rj) have been given earlier in this report. 

This correction technique is simple to use. Thus, it may be possible to use rela- 
tively large heat flux meters in obtaining acceptable local data. 

For example, for the Tabakoff and Clevenger (ref. 5) data, using equation (32) in 
conjunction with equation (30) gives a better estimate of the true value of the actual heat 
transfer coefficient at the stagnation point. The calculated value at the stagnation point 
is labeled "calculated peak" in figure 8. 

CONCLUDING REMARKS 

General formulas are derived for determining gage averaging errors of strip-type 
heat flux meters used in the measurement of one- dimensional heat flux distribution. In 
addition, a correction procedure is presented which allows a better estimate of the true 
value of the local heat flux. As an example of the technique, the formulas are applied to 
the cases of heat transfer to air slot jets impinging on flat and concave surfaces. For 
many practical problems, it is possible to use a relatively large gage to obtain accepta- 
ble local heat flux measurements, provided that the gage is small enough to detect any 
peaks which might occur in the heat flux distribution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 3, 1973, 

501-04. 
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Figure 1. - Relation between local and "average" heat flux 
for nonlinear heat flux. 




Figure 2. - Effect of gage size on measured error. 



Dimensionless distance from stagnation point, x/B 


(a) Impinging jet. (b) Wall jet. 

Figure 3. - Lateral variation of local heat transfer coefficient between a plate and 
an impinging two-dimensional air jet. Width of slot jet, B, 0.3175 centimeter 
(1/8 in.). (Data from ref. 4. ) 


y 




Figure 4. - Heat transfer for impinging 
hole jets on a concave surface. Hole 
diameter, D, 0. 635 centimeter 
(1/4 in.); spaci ng-to-diameter ratio, 
H/D, 1; height-to-diameter ratio, 
Z n /D, 8.8. (Data from ref. 5.) 
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Calculated peak 



Figure 8. - Heat transfer for row of holes impinging 
on concave surface. (Data from ref. 5.) 
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